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Abstract 

In order to investigate the genesis of powerful radio jet, we have mapped the central 10 pc region of the 
nearby radio galaxy NGC 6251 with a 0.2 pc resolution using Very Long Baseline Interferometer (VLBI) at 
two radio frequencies, 5 GHz and 15 GHz, we have found the sub-parsec-scale counterjet for the first time 
in this radio galaxy. This discovery allows us to investigate the jet acceleration based on the relativistic 
beaming model. 
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1. Introduction 

The powerful radio jets are observed in approximately 
ten percent of luminous AGN and their maximum lengths 
sometimes amount to ~ 1 Mpc (Bridle, Perley 1984). Al- 
though global morphological properties give us very im- 
portant information, very inner regions in the radio jets 
also provide hints for understanding the genesis of radio 
jets. The capability of Very Long Baseline Interferom- 
etry (VLBI) at radio frequency has been very useful in 
investigating such inner regions of relatively nearby radio 
galaxies (Zensus 1997). It has grown significantly since 
the launch of the space radio telescope HALCA (Highly 
Advanced Laboratory for Communications and Astron- 
omy); a key element of the VLBI Space Observatory Pro- 
gram (VSOP) system (Hirabayashi et al. 1998; this vol- 
ume). This has enabled us to investigate radio galaxies 
at the sub-milli arcsecond angular resolution, although 
the resolution depends on the radio frequency. 

Since it is known that powerful radio jets are often 
emanated from the central engine of AGN at a relativistic 
speed, the inner region of radio jets is just a very exciting 
place. Since a pair of radio lobes are generally observed 
in the radio galaxies, the jet is believed to be intrinsically 
ejected in two opposite directions. However, although the 
approaching jet is enhanced by the Doppler beaming, the 
receding one (the counterjet) is significantly extinguished 
if the jet velocity is highly relativistic and is observed 
from a smaller viewing angle toward the jet. Indeed, in 
many radio galaxies, the counterjet is hardly seen in the 
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vicinity of the radio core. Since a typical distance of the 
radio jets resolved by VLBI is ~ 1 pc inner region, this 
inner region is called as pc-scale jets (Zensus 1997). 

NGC 6251 is one of apparently brightest powerful ra- 
dio galaxies in the nearby universe and thus has been in- 
vestigated extensively using VLBI as well as VLA (Very 
Large Array) (Waggett et al. 1977; Cohen, Readhead 
1979; Perley et al. 1984; Jones et al. 1986; Jones, Wehrle 
1994). So far these measurements failed to detect the pc- 
scale counterjet. If we could detect it, it will be possible 
to give many important observational constraints on the 
jet geometry and then the physical process of jet accel- 
eration. If the jet acceleration would occur at a sub-pc 
region, we could detect the counterjet because the fad- 
ing due to the Doppler beaming is still insignificant. In 
order to find evidence for such sub-pc-scale counterjet in 
NGC 6251, we have performed new high-resolution VLBI 
observations using HALCA. We use a distance to NGC 
6251, 94.4 Mpc, which is determined with the use of a 
recession velocity to the galactic standard of rest, Vgsr 
= 7079 km 8"^ (de Vaucouleurs et al. 1991) and a Hub- 
ble constant Hq — 75 km Mpc~^. Note that 1 mas 
(milli arcsecond) corresponds to 0.48 pc at this distance. 

2. Observations 

NGC 6251 was observed at 5 GHz using VSOP on 
30 April 1998 (10.1 h on-source) and at 15 GHz using 
VLBA (Very Long Baseline Array) on 2 June 1998 (9.3 
h on-source) . The VSOP system used in our observation 
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Fig. 1.. The uv coverages at 5 GHz (a) and at 15 GHz (b). 



consists of HALCA, VLBA, and Effelsberg 100-m radio 
telescope near Bonn. Details of the observations are sum- 
marized in table 1. Since quality of the data taken at the 
Mauna Kea station of VLBA appeared to be poorer than 
those taken at the other stations, we flagged out the vis- 
ibilities with the Mauna Kea station. The band width 
of the VSOP receivers is 32 MHz, divided into two IF 
bands, each of which has 128 channels at 5 GHz. The 
bandwidth of the VLBA receivers is also 32 MHz but 
divided into four IF bands each of which has 16 chan- 
nels at 15 GHz. We used the source 1803-1-784 as a cal- 
ibrator at both frequencies. Correlation processes have 
been carried out using the VLBA correlater at Soccoro in 
New Mexico. The integration time at the correlater was 
3 s. We used the NRAO AIPS package for fringe fitting 
and a-priori amplitude calibration. Image restoration via 
self-calibrations were made using DIFMAP. The result- 
ing visibilities were coherently averaged for 30 s. The 
synthesized beam size of the two images are 0.73 x 0.34 
mas in PA = -9.2° at 5 GHz and 0.46 x 0.38 mas in PA 
= —13.4° at 15 GHz with uniform weighting. In order to 
perform beam-size-matched comparison between 5 GHz 
and 15 GHz, we restored the two images with a same 
spatial resolution of 0.50 x 0.50 mas. The uv coverages 
of the two observations are shown in figure 1. 

3. Results 

Our final maps at 5 GHz and 15 GHz are shown in 
figure 2a and 2b, respectively. The dynamic ranges are 
160 for the 5 GHz map and 1400 for the 15 GHz map. In 
both maps, several knots are found within 10 mas (4.8 
pc) from the each brightest peak. Analyzing the wiggling 
patterns in both frequencies, we obtain an opening angle 
of the jet, (/) ~ 0.50° ± 0.29°. It is remarkable that a 
counterjet-like feature can be seen in both maps although 
its spatial extension is only « 1 mas. The most prominent 
knot in the 5 GHz map is seen at w 2.0 mas (0.96 pc) 
away from the 5 GHz brightest peak. Its intensity is ~ 



26% of the 5 GHz core. This knot is found for the first 
time. It is also seen in the 15 GHz map. However, its 
angular distance from the 15 GHz brightest peak is ~ 2.3 
mas which is larger by 0.3 mas (0.14 pc) than that in the 
5 GHz map. 

Since any VLBI images have no information on the ab- 
solute position because of the self-calibration procedure, 
this difference makes it difficiilt to register the two im- 
ages precisely. Since it is often observed that radio cores 
are optically thick even at those frequencies because of 
higher plasma densities, it is imlikely that the observed 
brightest peak positions always correspond to the real 
core position; here we assume that the very center of the 
active nucleus, i.e., the central engine is enshrouded by 
the VLBI core. On the other hand, knots are thought to 
be optically thin and thus their brightness distributions 
show more realistic distributions of plasma. Therefore, 
it is better to register the two images using the spatial 
and morphological information of knots between the two 
maps. 

Here we apply two new methods for the registration. 
The first one is to analyze the correlation function be- 
tween intensity profiles of the radio jet at the two fre- 
quencies. The correlation function exhibits the maxi- 
mum correlation coefficient Cint — 0.98 at an offset of 
Aint — 0.30± 0.0025 mas, indicating that the peak at 
5 GHz is shifted by Aint to the jet direction from that 
at 15 GHz (figure 3a). The second method is to ana- 
lyze the correlation function between wiggle patterns of 
the radio jet at the two frequencies. This analysis gives 
an offset of A wig — 0.10 ± 0.24 mas with the maximum 
correlation coefficient of Cwig — 0.34 (figure 3b). Since 
an offset of Aint is in good agreement with that of A^ig 
within their errors, we adopt a A = 0.3 mas as the best 
offset. Using this offset, we register the two images and 
make a spectral index image of the radio jet of NGC 6251 
shown in figure 2c. Here the spectral index a is defined 
as Si, cc . It is shown that the core region has an 
inverted spectrum with a ~ 1.2 and thus it appears opti- 
cally thick. On the other hand, the spectral index of the 
jet ranges —1 and —0.5 and thus the jet is optically thin. 
The counterjet-like feature has a spectral index of ~ —0.5 
and thus appears to be optically thin. This ensures that 
it is the counterjet. Its intensity level is ^ lOfTims at 5 
GHz and ~ 40fTrms at 15 GHz, where (Jrms is the one 
sigma root-mean-square noise in each map. Therefore, 
we conclude that this is the real counterjet. Curiously, 
the morphology of the counterjet is slightly different be- 
tween the two images, i.e., two ridge lines can be seen 
only at 15 GHz. Although it is difficult to imderstand 
such a morphological difference, we confirmed that one 
of ridge lines at 15 GHz corresponds to the counterjet at 
5 GHz. 
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Table 1. Observations. 
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Peak Int. 


RMS noise 


Dynamic 




(GHz) 




(Jy/beam) 


(m Jy/beam) 


Range 


1998 Apr. 30 


5 


VLBA, HALCA*, ERt 


0.13 


0.82 


160 


1998 Jun. 02 


15 


VLBA 


0.34 


0.24 


1400 



'Highly Advanced Laboratory of Communication and Astronomy 
^Effelsberg 100-m telescope 
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Fig. 2.. The images of NGC 6251 at 5 GHz (a) and at 15 GHz 
(b). The spectral index map together with the 15 GHz contour 
image is also shown in (c). All the images rotated clockwise on 
the sky by 28.4°. 




Fig. 3.. Top: Intensity profiles (a) and wiggling patterns (b) at 5 
GHz (grey circle) and 15 GHz (filled circle). They are already 
registered. Bottom: Correlation coefficient as a function of 
offset. 



4. Discussion 

J^.l. The Absorption Model at the Core 

The inverted spectrum in the core implies that there is 
a strong absorption at 5 GHz. This may explain why the 
brightest position at 5 GHz does not coincide with that 
at 15 GHz. In general, such absorption in the radio jet 
is caused either by free-free absorption or by synchrotron 
self-absorption. Since it seems quite likely that a plasma 
halo (or ring) surrounds the central engine, the free-free 
absorption seems to be a dominant absorption mecha- 
nism, i.e., the flux density at 5 GHz is more significantly 
absorbed by the free- free absorption than the 15 GHz 
brightest peak flux. 

Free-free absorption models give a relation Si, = 
Soiy"" ex-p[—TQi^^^-^]. Accordingly, we obtain, 

_ - QQ / N -2.1 

0.0279 VgHzJ ' ^ ' 



where 015^ is the spectral index between 5 and 15 GHz 
and ao is the intrinsic spectral index of the jets. Using 
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2.4. 



1.2, 



ao 



-0.5, and ly — 5 GHz, we obtain Tobs 



A question arises as how the absorbing material is 
distributed around the central engine at a sub-pc-scale. 
Since we do not know the true situation, we will investi- 
gate the following three absorption models: (1) a spher- 
ical absorber model, (2) an optically-thick disk model, 
and (3) an optically-thin disk model. 

4.1.1. The spherical absorber model 

First, we investigate a spherical absorber model be- 
cause this model is the simplest one. We assume that 
the core is surrounded by a plasma sphere with a radius 
of a and the radio emission from the inner part of the 
jets (i.e., both the jet and the counterjet) suffer from the 
free-free absorption. We show this model schematically 
in figure 4. In this model, we adopt the 15 GHz bright- 
est peak as the core in which the central engine resides. 
The approaching jet escapes from the plasma sphere at 
a projected distance of a; = a x sin O^at and thus the op- 
tical depth becomes to be small (i.e., Tobs < 1) here. 
On the other hand, although the counterjet escapes the 
sphere at a; = — a x sin6'jot, it suffers the effect of free- 
free absorption until x — —a. It is also noted that the 
path length is longest at a; = — a x sin^jet . As shown 
in figure 4, we define the following projected distances; 



X 



jot 



X 



peak 



a x sinSjct and X, 



i-cjct 



Then we are 



counterjet 
(cjet) 




the core 

absorption 
region 



Fig. 4.. A schematic illustration of the spherical absorber model 
in NGC 6251. 
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In figure 5, we show the observed spectral index vari- 
ation along the radio jet. The opacity is assumed to 
be proportional to the path length through the plasma 
sphere. The region with an inverted spectrum is con- 
sidered to be optically thick for the free-free absorption 
{Toha > 1) while that with a normal spectrum is optically 
thin. We adopted that Tobs = 1 when the jets escape from 
the sphere. From equation (1), note that Tobs = 1 at 5 
GHz corresponds to — 0.2 and that Tobs at 15 GHz is 
always smaller than unity in our observations. Therefore, 
from figure 5, we estimate Xp^ak ~ 0.20 mas, Xjet ~ 0.25 
mas, and Xcjot ~ 0.75 mas. Then we obtain O-^et — 19°. 
In the previous VLBI measurements of the pc-scale jet 
of NGC 6251, Jones (1986) obtained 6'jct < 45° using 
the upper limit of the counterjet flux. Our estimate ap- 
pears consistent with their result. Although our simple 
model suggests A"jot = A"pcak, the observation results in 
ATjet/ATpoak ~ 1-3 7^ 1. This may be caused by plasma 
density fluctuation in the actual radio core of NGC 6251. 
If we adopt the observed value of X-^ct instead of A'poak, 
we obtain 6'jot — 15°. 




MilliARC SEC 



Fig. 5.. Spectral index variation along the jet near the core. 
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4-. 1.2. The optically-thick disk model 

Although the spherical absorber model is useful for 
estimating the jet viewing angle, it is too simple to de- 
termine the detailed environment of the radio core. Re- 
cently, it has been reported that a very narrow gap near 
the core of the radio jet in NGC 4261 is attributed to the 
absorption by a geometrically thin, nearly edge-on accre- 
tion disk (Jones et al. 2000). In this case, the central en- 
gine may be powered by gas accretion at a sub-Eddington 
accretion rate. It is expected that an accretion structure 
consists of an optically-thin, advection-dominated accre- 
tion flow (ADAF) in the central region and a geomet- 
rically thin accretion disk surrounding it (e.g., Narayan 
et al. 1998). This idea can be applicable to the case of 
NGC 6251. 

We consider that the central black hole is surrounded 
by such a tenuous plasma or coronal region which is re- 
lated to ADAF, and that a relatively dense outer disk 
causes a strong free-free absorption. In this model, the 
strongest absorption should be observed at the inner edge 
of the outer disk. Therefore, we identify the location of 
the maximum absorption at 5 GHz with the inner edge. 
If the jet viewing angle 0jct = 30° (see Sudou, Taniguchi 
2000), then we obtain the linear distance of the inner 
edge, Tpoak = Xpeak/ COS 6 « 3.3 X 10^'' cm. This model 
is shown schematically in figure 6. 

According to the usual understanding of accretion 
disks (Narayan et al. 1998), we regard the outer disk as 
the optically-thick disk (Shakura, Sunyaev 1973). Since 
the density and the temperature are low in this outer 
disk, the free-free absorption may be dominant over 
the electron 
(1973), the I 
sity n, and t 

H = 

n ~ 

and 

T = 

respectively, 
m is the bla 
and f is the 

are scaled fo. „ .j^. ^„ --^i,, ■■- — 

M/Mcrit (McHt = LEdd/0.1c2 ^ 1.5 X lO^^m g s~i), and 
r = r/rg (rg = 2GM/c^ ~ 3.0 x lO^^m cm). For NGC 
6251, we adopt m = 8 (Ferrarese, Ford 1999), m — 0.001 
(Birkinshaw, Worral 1994), and r — fpoak — 1-4 x 10^. 
Assuming avis = 0.1, we obtain, H = 4.0 x 10^"* cm, n = 
7.0 X 10^3 and T = 5.5 x 10^ K. Although H and 

T are reasonable compared with observations, n is too 
large as to given a very high opacity of r ~ 10^^ at 5 GHz. 
Therefore, this model fails to explain the observational 
properties of NGC 6251. 




Fig. 6.. A schematic illustration of the absorption disk model in 
NGC 6251. 



4-. 1.3. The optically-thin disk model 

Next, we attempt to explain Tobs adopting an optically- 
thin disk model in which the viscous heating is balanced 
by the free-free cooling. Adopting the cooling rate of 
1.2 X IQ'^'^mln'^T^/'^H erg s^^ cm-^ (Kato et al. 1998), 
we obtain, 

H = 5.8 X IQ^^a-^J^m^'^m^f'"^ cm, (7) 
n = 5.3 X lQ^^aJ.^m^''^m-^f-^^/^ cuT^ , an^) 
T = 1.0 X 10^'^a~^^rn^''^m^f~^/^ K. (9) 



4-. 1.4- Summary of the absorption models 

The simplest spherical absorption model provides the 
plausible explanation for the observed absorption proper- 
ties. On the other hand, the optically-thick and optically- 
thin disk models fail to explain the observations, al- 
though the presence of the disk structure would be nat- 
urally expected around the central engine. It will be 
necessary to seek for other types of outer accretion-disk 
models which are consistent with the observations. 



scattering. Following Shakura and Sunyaev 
lalf thickness of the disk H, the particle den- 
he electron temperature T are given by, 



4.0 X 10^°a 



-1/10^ 
vis 



,3/20^9/10-9/8 



2.1 X 10^'a 



m ' m ' r 

21„ -7/10^,11/20 -7/10 --15/8 



9.. 



10' 



1/5--3/4 ^ 



where avis is the viscosity parameter, and 
,ck hole mass, m is the mass accretion rate, 
radial distance. The latter three parameters 

m = rn — 



Adopting the same parameter values as those in section 
4.1.2, we obtain H ^ 9.0 x 10^^ cm, n = 4.7 x 10^ 
cm-3, and T = 1.4 x 10^ K for NGC 6251. The derived 
temperature of T ~ 10^ K is much higher than that of 
^ 10"* K assumed in the model of Jones et al. (2000). 
(3)i 

n view of the resulting high electron temperature, the 
eni~(i:,)sooling in this situation may be dominated by the inverse 
j-g-jCompton scattering rather than by free-free radiation. 

The model suitable for such a condition may be that of 
(6)Shapiro et al. (1976) which adopts that the radiation 
cooling is dominated by inverse Compton scattering. In 
fact, the expected free- free optical depth in our model is 
only ^ 10~®. Therefore, we conclude that this model is 
unsuitable for the absorption mechanism of NGC 6251. 
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4-2. Velocity of the Jet 

Adopting the so-called Doppler beaming model (e.g., 
Ghisellini et al. 1993), we estimate the jet velocity using 
the observed jet-counterjet intensity ratio; i.e., if the jet 
of NGC 6251 is intrinsically two-directional symmetric, 
the observed intensity asymmetry between the jet and 
counterjet can be attributed to the Doppler beaming. 
The jet-counterjet intensity ratio R is related to both the 
intrinsic jet velocity Vj^t and the viewing angle 0jet toward 
the jet by the formula 



R = 



1 + /?jet cos ^jet 



1 -/3jetC0s6ljet 



3-ao 



(10) 



where /3jci= Vjci/c. Using the intensity profiles of both 5 
GHz and 15 GHz along the jet direction, we have mea- 
sured R and plot them in figure 7. Note that the core 
component (|.x| < 0.5 mas) is not used in this analy- 
sis. Although the estimated value of R at 1.05 mas from 
the core is the same at both frequencies, that within 0.8 
mas from the core is higher at 5 GHz than at 15 GHz. 
Since this is probably due to severer absorption at 5 GHz 
than at 15 GHz, the value of i? at 15 GHz appears to be 
more reliable. It is shown that R increases with increas- 
ing projected distance from the core. If this increase is 
caused mainly by the Doppler beaming, it is suggested 
that /?jct, and/or 6jci increase with increasing distance. 
Since it is unlikely that 6'jet varies significantly at such 
a small scale, it seems reasonable to conclude that the 
increase of R is attributed to the increase in jSjct- This 
can be interpreted as possible evidence for the accelera- 
tion at the sub-pc-scale radio jet. Using equation (9) and 
assuming 0jct= 30 °, we show that the jet is accelerated 
from /3jet ~ 0.13 at 0.55 mas to /3jet ~ 0.42 at 1.05 mas. 
The results are summarized in table 2. These angular 
distances correspond to the linear distances of w 0.30 pc 
and « 0.58 pc, respectively. 

If the jet acceleration is effectively made at sub-pc- 
scalc region, the jet velocity becomes relativistic at radial 
distance of ~ 1 pc. If the sub-pc-scale acceleration is a 
general property in most radio galaxies, it seems difficult 
to detect pc-scale counterjcts. It will be very important 
to observe sub-pc-scale jets at higher frequency because 
the effect of absorption is negligibly small. 

Since the kinetic energy of a particle with a mass 
of TO at a relativistic velocity of v is given as E = 
rwyv^ = m^yP^cP where 7 = (1 — /?^)~^/^, the kinetic 
energy gain from /9 = 0.1to/3 = 0.4 per unit mass is 
AE/m ~ 3 X 10^^ J kg~^ This corresponds to the ki- 
netic energy of ~ 200 keV for an electron and to that of 
~ 300 MeV for a proton. Therefore, it is suggested that 
the jet acceleration region found in NGC 6251 is related 
to the 7-ray emitting region. This may explain why 7-ray 
loud AGNs tend to be associated with radio-loud AGNs 
if such acceleration occurs commonly in radio galaxies 
(Shrader et al. 1995). 
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Fig. 7.. The jet-counterjet ratio R as a. function of the projected 
distance from the core. The filled triangle and open circle in- 
dicate obtained R at 5 GHz and 15 GHz, respectively. 



5. Conclusions 

We discussed high-resolution images of NGC 6251 at 
5 GHz and 15 GHz. From the analysises of the inten- 
sity profile and the wiggling pattern, we register the two 
images and obtain the spectral index image. It shows 
the counterjet which has an optically-thin spectrum. In 
order to explain the position difference between the core 
and the absorption peak, we adopt the simple spheri- 
cal absorption model, and show the possibility that the 
jet viewing angle can be obtained from only the spec- 
tral index distribution. Although we also attempt to 
adopt the free- free absorption by an accretion disk model, 
both optically thick and thin disk models based on the 
a-viscosity cannot explain the observed opacity. 

The jet-counterjet intensity ratio analysis shows that 
the radio jet of NGC 6251 is accelerating from 0.13 c at 
the linear radial distance of r w 0.30 pc to 0.42 c at r » 
0.58 pc. This provides the first direct evidence for the 
sub-pc-scale acceleration of the radio jet in NGC 6251. 
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Table 2. The jet/counterjet intensity ratio. 



Distance (mas) 




Rl 


4t 


0.55 


6.8±0.4 


2.2±0.015 


0.13±0.0014 


0.80 


9.0±1.1 


6.3±0.23 


0.30±0.0053 


1.05 


13.5±3.0 


13.5±1.6 


0.42±0.017 



*R at 5 GHz 
at 15 GHz 
* Assuming ^jet = 30° 
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